Continuous quantum measurement: inelastic tunnelling and lack of current oscillations 
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We study the dynamics of a charge qubit, consisting of a single electron in a double well potential 
coupled to a point-contact (PC) electrometer, using the quantum trajectories formalism. Contrary 
to previous predictions, we show formally that, in the sub-Zeno limit, coherent oscillations in the 
detector output are suppressed, and the dynamics is dominated by inelastic processes in the PC. 
Furthermore, these reduce the detector efhciency and induce relaxation even when the source-drain 
bicis is zero. This is of practical significance since it means the detector will act as a source of 
decoherence. Finally, we show that the sub-Zeno dynamics is divided into two regimes: low- and 
high-bias in which the PC current power spectra show markedly different behaviour. 

PACS numbers: 73.63.Kv 85.35.Be, 03.65.Ta, 03.67.Lx 



Single shot quantum measurement of mcsoscopic sys- 
tems is recognised as an important goal. Fundamen- 
tally, it will allow us to make time-resolved observa- 
tions of quantum mechanical effects in such systems, and 
practically it will be a necessary component in the con- 
struction of solid-state quantum information processors 
(QIP). There are numerous proposals for implementing 
QIP in soHd state systems, using doped silicon 0, elec- 
trostatically defined quantum dots [21 and superconduct- 
ing boxes [3j. In these proposals, the output of the QIP 
is determined by measuring the position of a single elec- 
tron or Cooper pair. Single electrons hopping onto single 
quantum dots have been observed on a microsecond time 
scale using single-electron transistors (SET) Ensem- 




FIG. 1: Schematic of the qubit and PC showing lead energy 
bands. Electrons tunnelling from the source (S) to the drain 
(D) may do so elastically or inelastically, depicted by arrows. 
Different transitions induce different jumps, Pi, on the qubit. 



ble measurements of a double well system (qubit) have 
been demonstrated in superconducting devices [5j- So 
far, single shot qubit measurements remain elusive. 

It is therefore important to consider the measurement 
of single electron qubits by sensitive electrometers. Two 
possible electrometers have been discussed to date: SETs 
e.g. H) and point contacts (PCs) 0, S IU|,0 ElEl 



liall PCs are sensitive charge detectors|l5[ ll( 

and are the focus of this Letter. Figure ^ illustrates the 
physical system we consider here, with a PC (two Fermi 
seas separated by a tunnel barrier) in close proximity to 
one of the dots. 

Three energy scales are relevant: the qubit level split- 
ting, (j), the PC bias voltage, eV, and the measurement 
induced dephasing rate, F^ = 2{y/Ti — -^77)^/6, due to 
the distinct currents, Ii^r-, through the PC when the qubit 
electron is held in the left (1) or right (r) well. We ignore 
environmental decoherence [l^ . These energies define 
three distinct measurement regimes: (1) low-bias regime, 
Td/2 <t: eV < (j), (2) high-bias regime, Td/2 < < eV^ 
and (3) quantum Zeno Hmit, (f> <C Fd/2 ^ eV. 

In the quantum Zeno limit, frequent weak measure- 
ments localise the qubit, suppressing its dynamics |^linl|. 
To date, little attention has been paid to the low bias 
regime. Recently, the asymmetric power spectrum was 



calculated for arbitrary eV [l9j. This may correspond 
to the emission and absorption power spectra of quanta 
in the PC, as discussed in [i^- It has been predicted 
that in the limit, Td/2 ^ cp, referred to as the sub-Zeno 
limit, coherent oscillations due to the quantum dynamics 
of the qubit will be observed in the current through the 
PC nil [H [13 • It was also claimed that the PC is an 
efficient detector, i.e. no information about the qubit is 
lost by the detector. We present a detailed analysis of the 
sub-Zeno limit, showing the sharp difference in behaviour 
between the low- and high- bias regimes, and that previ- 
ously predicted coherent oscillations are suppressed. The 
analysis also formally yields the boundary at which ap- 
proximations leading to the Zeno effect are reasonable. 



For a low transparency PC, the Hamiltonians, in units 
where h ~ 1, for the qubit, leads, and their interaction 
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are given by H 0, i, @, E EEl 

i7,ys = (-Aa, -ea.)/2 = -(/.(Ti^)/2, (1) 

+ Xfc,gCT^)a^ ,^05,^ + H.c. (2) 

H\cads ~ Hs + Hd = ^ '^fc(agj.a5_fe + j,aD,fc) (3) 

k 

where a, = |0(r| + \r){l\, a, = - |r)(r|, ai"^^ = 

|5)(.g| - |e)(e|, = tan-i(#), = VA^T^, I5) = 
cos(6l/2)|0 + sin(6'/2)|r) and |e) = - sin(6'/2)|Z) + 
cos(0/2)|r). We will adopt the convention that Xk,q < 
so that the left well is nearest the PC. Furthermore we 
assume the PC is weakly responding so |xoo| |Tbo|- 

Unconditional master equation. The von Neumann 
equation for the density matrix, i?, of the bath and 
qubit is R(t) = — z[i7Tot,^]- To derive the master equa- 
tion (ME) for the reduced density matrix, p, of the 
qubit we transform to an interaction picture with re- 
spect to the free Hamiltonian Hq = Hgys + -ff leads, 
Hi{t) = e'"°^H^ei^e~'"°K We formally integrate the 
von Neumann equation, then trace over lead modes |2ll |. 

Mt) = TTsM-mit), R{0)h I dt'[Hi{t), [Hi{t'),R{t')]]}, 

Jo 

(4) 

where Hi{t) = Y.k.q^k,q{t)a^jj^^as,k + H.c, and 
Sk,q{t) = e-'("''-"')*(Tfc,, + Xk^q^rS))- We write 
Sk,q{t) as a discrete Fourier decomposition Sk,q{t) ~ 
g-it(<^fc-c^,)(g-it0p^ git^p^ p^)^ vjhcvc Pi = pI = 

-Xoosin(0)|5)(e|, P3 = (Too + Xoo cos(0)(Ti^'). We have 
assumed that Tk.q = Too and Xk,q ~ Xoo are constant. 

The form of Sk,q{t) indicates that there are three pos- 
sible jump processes, indicated in Fig.^ P3 is associated 
with elastic tunnelling of electrons through the PC. Pi 
(P2) is associated with inelastic excitation (relaxation) 
of electrons tunnelling through the PC with an energy 
transfer (jj. This energy is provided by the qubit which 
relaxes (excites) in response. Inelastic transitions in sim- 
ilar systems have been described in , which calculated 
the current power spectrum through an open double well 
system due to shot noise through a nearby PC. 

We make standard assumptions that result in a 
time-independent Markovian ME. Firstly we assume 
that the leads are always near thermal equilibrium 
so Trs,D{alSi?(i)} = and Trs,D{4fca„,feP(i)} = 
5i.jfi{'^k)pit), where i,j G {S', Z?}, fi is the Fermi distri- 
bution for lead i and p is the qubit density matrix. Sec- 
ondly, if the lead correlation time is much less than other 
time scales, then the lower limit on the integral in Eq. @ 
can be set to —00 We substitute J duj^gsio) 

where gs{D) are the lead density of states, which are as- 
sumed constant over the relevant energies. 

Finally, we make the rotating wave approximation 
(RWA), setting to zero terms in the ME proportional to 



g±i0t^ The RWA is apphcablc when > z/^eV (defining 
the sub-Zeno limit), where v — ^/^ngsgoXoo, which is 
small in the low transparency, weakly responding limit. 
We also make the approximation p{t') — > p{t) in Eq. Q), 
which is valid when 1, to arrive at the interaction 

picture ME 

pi(t) = 2TTgsgD (T^[VWvT^Pn]piit) 

n 

+ VWe{-eV -u;n)P,\]piit)) = Cipiit), (5) 

where eV = /is — po is the bias applied across the PC, 
V[B]p = J[B]p - A[B]p, J[B]p = BpB\ A[B]p = 
\[B^Bp + pB^B), e{x) = (x + \x\)/2 , uJi = -u^ = 0, 

= and fii is the chemical potential of lead i. The first 
term in Eq. 10) arises from forward tunnelling electrons 
(S-D), whilst the second is due to reverse processes. We 
note that the dynamics implicit in this ME differs from 
those presented previously [a S H 0, ^| , since it includes 
inelastic jump processes. 

In the Zeno limit, ^ ly^eV, the RWA breaks down. 
Instead an alternate approximation is accurate, e**^* 1, 
or cquivalcntly ffzit) — s- ^^(O). In this limit, Sk,q{t) — 
g-i(wfc-w,)t(^P^^ -I- XooO'z(O)) has only one Fourier compo- 
nent, and we regain the ME of 0, 11,13 • 

Conditional ME. By unravelling the ME into evo- 
lution conditioned on measurement results, we can cal- 
culate the observed current correlation function and the 
power spectrum. From the unravelling it is also pos- 
sible to simulate quantum state trajectories consistent 
with single-shot experimental realisations of measured 
currents. There is no unique unravelling for Eq. one 
may add an arbitrary constant to each of the jump oper- 
ators in Eq. ©, to produce new (rescaled) operators 

ci^iy^eV + ^smid)\g)(e\+ji, (6) 
C2 = J^V|el^-0|sin(0)|e)(.g|+72, (7) 
C3 = i^VeV cos{e)ai''^ + Ty/eV, (8) 

which along with a modification to the system Hamilto- 
nian, provides an identical unconditional ME. Note that 
T = ^2TTgsgDToo. 

The unravelling that most accurately represents a 
given measurement process must be determined from 
physical considerations. We discuss in detail only the 
unravelling in the high bias regime, in which all jumps 
correspond to tunnelling from source (S) to drain (D), 
as reverse processes are Pauli blocked. Thus, when mea- 
suring currents through the leads, none of the processes 
are distinguished, since they all just contribute to a cur- 
rent. As a result, when a jump occurs, the resultant 
state of the qubit is a probabilistic mixture pic{t + dt) = 
X)n '^['^"]^'c(i)c?i- Between jumps, it evolves smoothly 
according to pocit) = -i[Hsys, poc{t)] - I]„ .4[c„]poc(i)- 

The constants associated with the inelastic processes, 
7i_2, are eliminated using energy conservation arguments. 
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Imagine we prepare the qubit in the state \g). Then sup- 
pose, using a sensitive bolometer measuring the change 
in energy in the leads, we determined that an inelastic 
jump had followed. We would conclude that a lead elec- 
tron tunnelling from S to D lost an amount of energy 
—0^2 = <p- Simultaneously, the qubit evolves discontin- 
uously, \g) C2|.g). If 72 = then this state is |e), so 
the qubit gains an energy 0, and energy is conserved, 
as required. Otherwise, if 72 7^ then it is straightfor- 
ward to show that energy is not conserved on average. 
A similar argument gives 71 = 0, and these arguments 
also hold in the low bias regime. The same reasoning 
applies when considering current measurements (rather 
than bolometric), resulting in 71 = 72 = 0. 

In order to determine the final constant in the unrav- 
elling, T, we relate the current through the PC to the 
jump processes according to the relation 

I{t)dt = e Tr {pi,(t + dt)] = e ^ Tr {J[cMt)} dt. 

n 

(9) 

In the absence of tunnelling {6 = 0), if the qubit is lo- 
calised in the energy eigenstate \l) (|r)) we expect a cur- 
rent Ii = eT^ [Ir = eT^) to flow through the PC. Solving 
the two equations Y.n Tr {c,i|0(%U = eV{T+vY = Tf 
and Y.n Tr {c„|r) (r|ct } = eV{T - vf = yields T = 
{Tr + Ti)/2VeV and v = [Ti - Tr)!^^^. The average 
and difference in the currents arc /loc — {Ir + Ii)/'^ = 
e(T2 + i/2)el/ and SIioc ^ U - h ^ ~4eiyTeV. 

For general values of 9 the mean current does not 
change, I {9) = I\oc, but the difference in current between 
the two localised states is SI{6) ~ (JJioc cos^(0), to low- 
est order in v. Notably, the variation in current between 
localised states vanishes when 9 = it/2. This demon- 
strates that when the energy eigenstates are completely 
delocalised, the PC measurement does not localise the 
qubit. This is evident from the form of the jump oper- 
ators, since when 9 = tt/2, C3 does not affect the qubit 
whilst C1.2 only induce transitions between the symmetric 
{\g)) and antisymmetric (|e)) states. 

Results. In the high bias regime, the ME is given by 
Eq. Q with Q{—eV — lj„) = 0. Solving the high-bias 
ME in the interaction picture gives 

p,e(i)=e-I^^(l+™^^(''))*/V,e(0), 

P.. W = ^ - ^(^)*(^ - P,,(0))aO) 

where we have defined Td ~ 2i^^eV. 

Our analysis is valid in the low bias regime, eV < 0, 
also. The unconditional ME is given by Eq. lO with 
e{-eV + 0) > 0, so it includes the term 'D[cl]p{t). The 
solution to the unconditional ME in the low bias regime 
is 

p,e(t)=e-(r^-^^W+^^^""'(«))V.e(0), 

p,,(t) = 1-6-2"^*^^"^ W*(l-p,,(0)), (11) 
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FIG. 2: Equilibrium ground state occupation probability 
for a qubit near a PC (solid), and for a thermalised qubit, 
Pgg°™(oo), at a temperature T = eV/2 (dashed). 



Note that a.t (f> ~ eV, Eqs. and ifTT)) agree. 

The unconditional, steady-state probability of the 
qubit to be found in its ground state is given by pgg (00) = 
{g\p{oo)\g). From Eqs. ((Tn|) and |[TT|) . we see that pgg(oo) 
depends on the regime in which the PC is operating 

In contrast, if the qubit were in thermal equilibrium with 
a heat bath at temperature T, then the ground state oc- 
cupation would be p^^™((X)) = e'^/^/(l e"^/^). Fig- 
ure 121 shows these two ground state probabilities when 
T ~ eV 12. There is an evident analogy between the PC 
bias voltage and an external heat bath. We note that the 
leads are each nominally at zero temperature. 

This correspondence between the PC bias and a tem- 
perature agrees with the predictions of [lH which shows 
that a PC induces effective thermal fluctuations in the 
position of a nearby harmonic oscillator. Even in the 
limit eV 0, the PC still induces relaxation, and 
Eq. (|11|) shows that in this regime, the energy relaxation 
time is = 2v^ (j) sin^ {9) , and the dephasing time is 
= TdCos^{9) + iy^(l)sin^{9). When eV" = or 6* = 7r/2, 
Td = 2Tr, indicating that the dephasing is solely due to 
energy relaxation, much like the decay of an atom in a 
vacuum. This is practically significant, as it means the 
PC cannot be turned off merely by making eV = 0. 

Another important time scale is the measurement time, 
r„i. Following 0, we calculate the initial rate at which 
the z-componcnt of the Bloch vector increases, starting 
from the symmetric state |^/'(0)) = {\l) + \r))/\/2. This is 
given by EiTiiai"^ dp,{t)y]/2 = r^^dt. We find r^^ = 
rrfcos^(^) -I- 0{v^(j)). Also, Tin > Td, with equality only 
when 9 = 0, indicating that the detector is inefficient 
unless the qubit energy eigenstates are localised. 

We now calculate current power spectra, using the two- 
time correlation function G(r) = E[I{t+T)I{t)]- E[I{t+ 
T)]E[I{t)], where E[...] is the classical expectation. Using 
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FIG. 3: Current power spectra, S(ijj) = ^^jj^{S(ijj) — 

S{oo)). For a given v the plot is valid in the region eV/ff) 
l/u^. Note the sharp change in the power spectrum at 
= ey. 



Figure 13 shows power spectra for different eV/cf). These 
power spectra are notably different from those predicted 
elsewhere in the sub-Zeno limit 0, 0, 0, 0|, with the 
absence of coherent oscillations at oj = ±0. 

In conclusion, we have shown that, in the sub-Zeno 
limit, inelastic tunnelling processes through a PC are 
important, and coherent oscillations are absent in the 
detector output, contrary to previous claims. Projective 
measurements in the localised basis are still possible as 
long as the energy eigenstates themselves are localised. If 
the eigenstates are not localised, inelastic jumps reduce 
the detector efficiency. These inelastic jumps also gener- 
ate pseudo-thermal fluctuations in the qubit. This is true 
even when eV = 0, leading to dccoherence of the qubit 
when the detector is nominally off. This is of practical 
significance in a QIP, since it means the detector will act 
as a source of decoherence. 



Eq. Q this can be written as |g 

G{T)^e\Tv{Y,J[cn]e''-'^^JK,]p^} - Tr{^ J7[c„]p„,}2) 
+ e2Tr{^^[c„]poo}5(T), (13) 

n 

where Cs is the Schrodinger picture version of the evo- 
lution operator in Eq. The power spectrum is then 
S{lo) = 2 dr G{T)e'"^^. Since /(t) is the real- valued, 
continuously-measured current through the PC, G(r) is 
manifestly symmetric. To compute the steady state cor- 
relation function, we take the limit t — > oo, so p(t) —^ poo. 
Keeping only the lowest order terms in a series expansion 
in (5/ioc, then taking the Fourier transform gives 



S'hb {uj) = Sq-c SIioc y cos{e) 



where So = 2e/ioc. 

In the low bias regime, two distinguishable jump pro- 
cesses occur: S-D electron tunnelling (elastic and inelas- 
tic), and D-S (inelastic). Therefore the current is re- 
lated to the number of jumps in the time interval by 
idt = e{dN+{t) - dN-{t)), where dN+^'^t) counts the 
number of S-D (D-S) tunnelling events in the time in- 
terval {t, t + dt). We find the steady-state, low bias cur- 
rent correlation function is just due to clastic tunnelling 
through the PC, 

5ib(co) = 5o - e «ioc cos(0) + 0(«iL). 
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